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Abstract 

We prove that any isotropic positive definite function on the sphere can be 
written as the spherical self-convolution of an isotropic real- valued function. It is 
known that isotropic positive definite functions on d-dimcnsional Euclidean space 
admit a continuous derivative of order [(d — l)/2]. We show that the same holds 
true for isotropic positive definite functions on spheres and prove that this result 
is optimal for all odd dimensions. 



1 Introduction 

For an integer d £ N we denote the d-dimensional unit sphere by S d = {x E R rf+1 
1}, where ||-|| denotes the Euclidean norm on R d+1 . A function f : S d x S d 
positive definite if 



is 



i=l j=l 

for all sets of points Ui, . . . , tt n £ § d and coefficients c±, . . . , c n £ R. The function / is 
isotropic if there exists a function / : [0, ir) — > R that fulfils 

f(u, v) = f(9(u, v)) for all u, v G S d , (2) 

where the geodesic distance on S d is given by 9 : §> d x § rf — > R, 0(u, f ) = arccos((u, u)). 
Here, (•, •) denotes the standard scalar product on IR^" 1 " 1 . 

Isotropic positive definite functions on spheres occur in statistics as correlation 
functions of homogeneous random fields on spheres or of star-shaped random parti- 
cles. They also have applications in approximation theory where they are used as 
radial basis functions for interpolating scattered data on spherical domains. Recent 
applications in spatial statistics can be found in Banerjee (2005); Huang et al. (2011); 
Hansen et al. (2011); application examples in approximation theory are given in the 
works of Xu and Cheney (1992); Fasshauer and Schumaker (1998); Cavoretto and Rossi 
(2010). 
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The class consists of all continuous functions tp : [0, tt] — > M. with ip(0) = 
1, such that the isotropic function ijj (6 (•,•)) is positive definite. The classes are 
nonincreasing in d, 

oo 
d=l 

with the inclusions being strict. 

We define the spherical convolution of two functions /, g : S d x E> d — > R as 

(f®9)( u i v ) = / f(u,w)g(w,v)dw, for all u, v G S d , 

where the integration is with respect to the d-dimensional Hausdorff measure on S d . 
The total measure of S d is denoted by a d = 2vr( d+1 )/ 2 /T((d + l)/2). It is easy to see 
that the spherical self-convolution of any symmetric L 2 -function / on § d x S rf is positive 
definite. 

Spherical convolution has been used by Wood (1995); Schreiner (1997); 
Estrade and Istas (2010); Hansen et al. (2011) as a tool to construct spherical posi- 
tive definite functions. It is natural to ask the reverse question: Which functions can 
be obtained through this construction principle? We can give the following general 
positive answer, which we prove in Section 3. 

Theorem 1.1. Any ip G has a spherical convolution root, which can be taken to be 
real-valued and isotropic. 

The techniques used to show the convolution representation theorem have lead to 
the solution of a further interesting problem concerning positive definite functions on 
spheres. 

A positive definite function / on the Euclidean space M. d is defined analogously 
to (1). The function / is called radial, if f(x,y) = f(\\x — y\\) for some function 
/ : [0, oo) — > M. Schoenberg (1938) showed that radial positive definite functions on M. d 
have a continuous derivative of order [(d— l)/2], where [c] denotes the greatest integer 
less or equal to c. The following theorem, which will be shown in Section 4.1 confirms 
the conjecture of Gneiting (2012) that the same holds true on spheres. 

Theorem 1.2. The functions in the class ^ admit a continuous derivative of order 
[(d — l)/2] on the open interval (0,7r). 

The derivatives at the point -d = can be infinite or can take finite values. We 
believe that the same holds true at $ = it. However, we are currently not able to 
provide simple examples for the latter claim. The powered exponential family 

^(0)=exp(- (*) a ), tfe[0,7r] 

with parameters c > and a G (0, 1] belongs to ^c^; see Gneiting (2012). For a < 1 
the first derivative at zero is — oo, whereas for a = 1 it takes the value — 1/c. The sine 
power function 



= 1 - 




■& G [0, vr] 
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of Soubeyrand et al. (2008) is a member of for a £ [0,2]. For a € (0,1), the 
first derivative at zero is — oo; for a = 1, we obtain ^'(O) = —1/2. If a € (1,2], the 
derivative at zero is zero. 

In the Euclidean case it is known that Theorem 1.2 is the best possible; see Gneiting 
(1999). Hence, there are radial positive definite functions on M. d whose derivative 
of order [(d — l)/2] + 1 is not continuous. The optimality of Theorem 1.2 for d = 
1,3,5,7 follows from the results of Beatson et al. (2011). In section 4.2 we introduce 
a turning bands operator for isotropic positive definite functions on spheres to show 
the optimality of Theorem 1.2 for all odd dimensions. In even dimensions it remains 
an open problem. However, once the optimality can be shown for d = 2, the turning 
bands operator immediately also yields the assertion in all even dimensions as well. 

The convolution representation result, Theorem 1.1, also has consequences that are 
of interest in statistical applications. Firstly, it shows, that any isotropic covariance 
function on the sphere can be obtained by the Levy based approach to modelling star- 
shaped random particles introduced by Hansen et al. (2011). Secondly, the proof of 
Theorem 1.1 reveals a way to resolve the identifyability issues associated with these 
models. It is possible to distinguish one specific convolution root amongst all possible 
convolution roots of a given covariance function. This is the basis of the inference 
procedure described in Ziegel (2012). 

2 Convolution of isotropic functions on spheres 

Let L 2 (E> d x § d ) be the space of square-integrable functions on x S d with the Hausdorff 
measure. By (■, -) L 2 and ||-||i2 we denote the scalar product and the norm of the Hilbert 
space L 2 (E> d x E> d ), respectively. We consider the closed subspace L 2 dX C L 2 (§> d x E> d ) 
of functions that are isotropic as defined at (2). For / € L dX it holds for all d + 1- 
dimensional orthogonal matrices R that 



f(Ru,Rv) = f(9(Ru,Rv)) = f(0(u,v)) 



f(u,v), for all u,v € S d . 



This property characterises the functions in l? d 



d,x 



Proposition 2.1. The convolution f ® g of f,g £ L^j- is in l? d 



d,X 



and 



||/®0lUa < a d sup \{f ®g){u,v)\ <\\f\\ L 2\\g\\ L 2. 



(3) 



u,v€§ d 



The convolution is bilinear, commutative and 



\\f®g\\h = (f®f>9®g)L*- 



(4) 
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Proof. It is easy to check that / © g is isotropic. Furthermore, by Holder's inequality, 

\(f®g)(u,v)\< [ \f(0(u,w))g(0(w,v))\dw 



1/2 f r 1 1/2 

< i I f(6(u,w)) 2 dw\ / g(6(w,v)) 2 dw 



d 



( = } /J- f [ f(e{u,w)) 2 dwdu\ l {—If g(6{w,v)) 2 dwdv\ ' 
{o-dJs d J§< i J [o-dJs d Js d J 

= — ll/H^IMI^ 

for u, v G S d . The equality at (*) holds true, because the integrals on the left hand 
side do not depend on u, v, respectively. Therefore, we obtain (3), and, in particular, 
/ © g G L 2 (S d x S d ). Bilinearity and commutativity are clear, and equation (4) is an 
application of Fubini's theorem. □ 

Schoenberg (1942) characterised the functions of the classes ^?d using Gegenbauer 
(or ultraspherical) polynomials. Let A > 0. The Gegenbauer polynomials for n G No 
are defined by the expansion 

1 oo 

— _ -^ = yVc£(cos0), for0€[O,7r]; 

(1 + r 2 - 2rcostf) A ^ 

v ; n=0 

see DLMF, 18.12.4. We will repeatedly use that 

A _ r(n + 2A) 
Cn{1) ~ n\T(2X) ■ (5) 

If A = we set C°(cos$) = cos(n#) for i? G [0, 7r] as in Schoenberg (1942). We need 
the following important property of the Gegenbauer polynomials with X = (d — l)/2, 
which follows from Xu (2005, Theorem 3.7). For d > 2, k,n G No and u, v G S d , we 
have 

Ct^d^w^f-^dw^^dw = S k>n a d - d ~ l c!t X)l \M), (6) 
d in + d — 1 

where 8k,n denotes the Kronecker delta. If A = it holds that 

C%((u,w))C°((w,v))dw = S k>n 7rC°((u,v)) 

for n G N , k G N, u, v G S d , and J gl Cg«u,«;))Cg((«7,t;»dti; = 2vr. 

Proposition 2.2. Let d > 2. The family = {Ed >n } n ^ , where Ed^ n '■= Cd,n x 
C { n d - 1]/ \(;-))eL 2 dI With 



2n + d - 1 



(d-l)C^- 1)/2 (l) 
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is an orthonormal basis of L? dX . Furthermore, for k,n € No, 

Ed,k © Ed, n = 0~k,nC~d,nEd,n, 

where 

C-d,n ~ 



d-1 



{2n + d-l)C [ n~ 1),2 {l) 

Proof. By (6) 

d-1 
l 2n + d - 1 



C^ l)/ \{uM)C^~ l),2 {{n,v))dudv = 5 Kn a d <2n ■ - - / C^- 1 )/ 2 ^,^ 

d-1 



X 2 



2n + d 



T cr i)/2 (i) 



hence is an orthonormal system. It is also a Hilbert space basis, because polynomials 
are dense in L 2 ([— 1, 1]). The second assertion is a direct consequence of (6). □ 

The following Proposition complements Proposition 2.2 and is not hard to prove. 
Proposition 2.3. Proposition 2.2 also holds for d = 1 with 

fl/(2vr), forn = 0, _fl, for n = 0, 

Cl '™ \y/2/(2ir), forn>l/ ° l,n [V2/2, forn>\. 

Propositions 2.2 and 2.3 imply that, for any function / £ Ljj, we have 

/ = (/; E dt n} L 2Ed :n , 

n£No 

where = means that the series on the right hand side converges unconditionally in L 2 
to the left hand side. We call the basis Cd the Gegenbauer basis of L dX . The coefficients 
(f,Ed t n)L 2 are termed the Gegenbauer coefficients of /. 

Proposition 2.4. For any f ^ L dX , n ^ No, we have 

f © Ed lU = C~d,n(f, Ed^n) L 2 Ed,n- 

Proof. For N G N we set /at = J2k=o(f> E d,k) L 2E d,k- Then /jy converges to / in L 2 . 
We obtain 

11/ © Ed >n - C~d,n{f,Ed, n )L 2 Ed,n\\L 2 

< 11/ © E d , n - /JV © ^d,n||L2 + \\f N ® E d (/) ^d,n)L2 J E (Jin || L 2. 

The last summand on the right hand side is zero by the definition of f?q and Proposition 
2.2. By Proposition 2.1 we obtain 

\\f ®E djn - f N ®E dtTl \\ L 2 = \\(f - f N ) ® E d , n \\ L 2 < \\f - fN\\ L 4E d ,n\\L 2 ->-0, 
as N — > oo. □ 
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Corollary 2.5. For any f G L 2 dX , n G No we have 

(f ® f,E d ,n)L2 =Cd,n{f, E d,n)h- 

Proof. We have 

{f®f,E din ) L 2 = (C djn ) _1 {/ © f,E dyn ® E djn ) L 2 = (c^n^Wf ® E d , n \\ 2 L2 

= (c d ,n)~ \\cd,n(f'E d ,n)L 2 E d)n \\ L 2 = C d ,n(f, E d)n ) L 2 , 

where we used Propositions 2.2 and 2.3, equation (4), and Proposition 2.4 in this 
order. □ 

The following theorem gives a necessary condition for the existence of convolution 
roots in L dX - In the interesting special case of nonnegative Gegenbauer coefficients 
this condition is also sufficient. 

Theorem 2.6. If a function f G L dX can be represented as f = g®g for some g G L dX 
then 

oo 

X^n^K/^nM <°°- ( 7 ) 
n=0 

// (7) ZioWs and (f,E djn ) L 2 > /or a// n G No, £/ien £/iere exists a g ^ L? dX such that 
f = 9®9- The coefficients of g in the Gegenbauer basis can be chosen to be nonnegative. 

Proof. The Hilbert space L dX is isometric to the space £ 2 (Werner, 2002, Corollary 
V.4.13). Therefore £] n6No a n E d ^ n G L? dX if and only if (a n )„ G N G or, equivalently, 
X]^Lo°n < 00 • Suppose now that / is given by / = g ® g for some g G -f^j- By 
Corollary 2.5 we have that 

{g,E d>n ) L 2 = ±(c d)n )- 1 / 2 |(/,^ ) „) i2 |V2 ) 

hence 

oo 

,u)lA < 00 • 

n=0 

For the reverse implication set g = J2neN ( d d,n)~ 1/2 (f, E^n} 1 /? E d>n . By assumption 
g G L dX and by Corollary 2.5 we have for any n G No, that 

(g © g,E d>n ) L 2 = c d , n (g,E djn ) L 2 = (f,E d)Tl ) L 2. 

With Parseval's equality (Werner, 2002, Theorem V.4.9) this yields the claim. □ 

We conclude this section with a proposition that shows that convolution products 
can be uniformly approximated with respect to the Gegenbauer basis C d . 



Proposition 2.7. If f G L dX is given by f = g ® g for some g G L d x> then for every 

permutation a : N — > N, 
converges uniformly to f . 



permutation a : N -> N, the sequence (/aOatgn with f N = £)fc=o(/> E d ^^} L 2E d ^^ 
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Proof. Let g N = Efe=o(f » E d,a{k)) tf E d,a(k)- By Corollary 2.5 and Proposition 2.4 we 
have 

N 

f ~ fN = 9® 9 Cd,o(k) (ff> E d,a(k))h E d,a(k) 

fc=0 
N 

= 9® 9- ^2(9 1 E d,a(k)}^9 © E d,a(k) = 9®g-g®9N = g®{g- 9n)- 

k=0 

Now, we can apply Proposition 2.1 to the last term and use the unconditional L 2 - 
convergence of to g in order to obtain the claim. □ 



3 Convolution roots 

Schoenberg's characterisation of the classes ^ is summarised in the following theorem; 
cf. Schoenberg (1942). 

Theorem 3.1 (Schoenberg). The class ^ consists of all functions of the form 

co fy(d— 1)/2/ q \ 

m = E^n Cn r(d farted), 

n=0 (1J 

with nonnegative coefficients b^ n > such that Yl^Lo ^d,n = 1- If d = 1, then 

6i = - / and b ln = - I cos(ra?Wtf)dtf, forn>\. (8) 

TT Jo TT 7 

If d>2, then for n £ N 

= 2n 2 tt" 1 rffif £ {ck'-v^m^^md*. 0) 

For a function ?/; £ vf^i we call the associated coefficients 6^ jn as given by (8) or (9), 
respectively, the d- dimensional Schoenberg coefficients of if). 

A function ip £ ^ is strictly positive definite if the inequality in (1) is strict for 
all systems of pairwise distinct points, unless all the coefficients are zero. Chen et al. 
(2003) show that ip £ ^ for d > 2 is strictly positive definite if and only if its 
Schoenberg coefficients bd, n are strictly positive for infinitely many even and infinitely 
many odd integers n. The corresponding result for is was derived by Menegatto 
(1994). Despite recent advances Sun (2005) there is no concise characterisation of the 
strictly positive definite functions in in terms of non-zero Schoenberg coefficients 
available. 

We prove the following result, which is slightly more detailed than Theorem 1.1. 
Theorem 3.2. For any ip £ ^ there exists a function g £ L\ x , such that 

ip(9(u, v)) = (g © g)(u, v), for all u, v £ E> d , 
and g has nonnegative Gegenbauer coefficients. 
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Proof. First, let d > 2, tp G The nonnegative Schoenberg coefficients of ip are 
connected to the Gegenbauer coefficients of ip(0(-, •)) via 

b d ,n = 2n 2 tt" 1 rff iy ^gr i)/2 (cos^)(sin^)^VW^ 

C^~ 1)/2 ((-u, v))i)(6(u, v))dudv 

S d xS d 

(r(^i)) 2 r(|)(d-i) ^ 

= r(d - l)2 4 - rf 7 r( rf + 1 )/ 2 (Cd ' n) ^'W - ' - )))^- 

The quotient in the previous line is positive and only depends on d. We denote it by 
ay. In particular, {E^ n , ip{0{--> -)))l 2 > f° r an n £ ^o- We have 

Ck d " 1)/2 ((v)) _ T - p 

hence 

00 

•)) = "rf^rf ^(^ n ,V(^(T)))L2^,n. 
n=0 

By Theorem 3.1 

00 00 00 

ra=0 n=0 ra=0 

hence Theorem 2.6 yields the claim. For d = 1 we have 



h/(2ir)(E lin ,i>(e(;-))) L 2, if n = 0, 
\A/2/(27r)(^, n ,^(fl(-, •)))£». if n>l, 



hence we can apply the same arguments as above. □ 

Remark. For a function tp G ^d+k C Vl/d for some k > 1, Theorem 3.2 yields spherical 
convolution roots (fa+fc G L 2 d+kX and ^ G L 2 ,^- with respect to the convolution in 
gd+fc anc j gd^ respectively. The associated functions g~d+k, 9d are both defined on [0, tt] 
and one would hope for a simple functional relationship between them, but it remains 
elusive thus far. However, on the level of Schoenberg coefficients, the functions gd+2 
and gd are easily put in relation using Gneiting (2012, Corollary 3). 

Let tp G ^d- The construction in the proofs of Theorems 2.6 and 3.2 shows that 
the class G d {ip) of all spherical convolution roots g G L d x of tp is given by all functions 
j £ L^j, whose Gegenbauer coefficients are given by 

K 1/2 <) n ,N . (io) 
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where (6d, n )neNo are the Schoenberg coefficients of ijj and (a n )n£N is a sequence with 
a n G {—1,1}; cf. Figure 1. In Theorem 3.2 we identify a unique convolution root 
by setting a n = 1 for all n 6 Nq. This choice resolves the identifyability issue when 
inferring the kernel of Levy based models for star-shaped random particles from their 
covariance or correlation structure as mentioned in Section 1. See also Hansen et al. 
(2011); Ziegel (2012). 

We conclude the section by using the convolution representation to calculate the 
Schoenberg coefficients of the function 

1 



t d : [0, tt] -> R,# M- —r-r 1 •) < r} © 1{0(-, •) < r}(0), 

where r G (0, vr/2], and is the normalising constant ensuring that i d (0) = 1. The 
convolution is taken in S d x S rf . It is a short calculation to show that v\{r) = 2r. For 
d > 2 the normalising constant is given by 

^(r) = a d _ 1 /' r (sin0) d - 1 d0. (11) 
J o 

The function £2 has been calculated explicitly by Tovchigrechko and Vakser (2001). 
Estrade and Istas (2010) provide a recursive formula for the functions L d , d > 2. 

Lemma 3.3. Lei r G (0, 7r/2]. The function !{#(•,•) < r} G L^ x has Gegenbauer 
coefficients {uJd,n}neN given, for n > 1, by 

Ud,n = c d ^ n o d o- d _i - 1 — (sin(r)) d C^ 1)/2 (cos(r)), for d>2, 
n[n + a — lj 

and ui\ )Tl = (2y/2/n) sin(nr). Finally, uj d fi = v d {r), where v d {r) is given at (11). 
Proof. Suppose first that d>2. We have 

(1{0( V ) ^r},^)^ = c d , re / / l{e{u,v)<r}C^' 2 {{u,v))dudv 
d JS d Js d 

(p7V \ 1 /*7T 

y (sin#) d ~ 1 d'&\ J t{#<r}Cb d - 1 V 2 (cos#)(sin'&) d - 1 d'& 

J cos(r) 

Using c d fl = cr^ 1 , the formula for n = follows. By DLMF, 18.9.20 we have for n > 1 



d_ 

d,r 



((l - *y/ 2 cS 1)/2 (*)) = - n( " d + -i" 1) (1 " * 2 ) ( "~ 2)/2 ^~ 1)/2 (*)> ( 12 ) 



which implies the lemma. The case d = 1 is a simple calculation. □ 

Using the relation between the Gegenbauer and the Schoenberg coefficients calcu- 
lated in the proof of Theorem 3.2 we obtain the following corollary. 
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o 
c\i 



D) O 



o 
o 



— ^ — ^ _ 



i 1 1 1 1 1 r~ 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 1: Different convolution roots g of L2( r ) for r = 1.2. The solid lines display 
the function z^(1.2) —1 / 2 l{i9 < 1.2} and its approximation by the first 32 Gegenbauer 
polynomials. The dashed line is the convolution root with nonnegative Gegenbauer 
coefficients. The dotted line represents the convolution root with a n = (— l) n , whereas 
the dash-dotted line has a n = (— ly n ' 2 \ with (o" n ) n , e N as in (10). 



Corollary 3.4. The function i& is in tD^. For d > 2 its Schoenberg coefficients are 
given by 

d '° aj T{d - l)24-<M<m)/2 ' 

and, for n > 1, 

/ r {d+l)/2, ,\ 2 



where 



Wr) = — — * -j sinrr". 



For d = 1, we /iaue &i o = r/(4-7r 3 ) and &i jn = v2 sin 2 (nr) /(rn 2 7r 2 ) /or n > 1. 

This example illustrates that the convolution root constructed in Theorem 3.2 may 
not be the most natural one. The Gegenbauer coefficients of i^(r) _1 / 2 l{#(-, •) < r} 
take both, positive and negative, signs; cf. Lemma 3.3. Hence, it is not the convolution 
root of id that results from the construction in Theorem 3.2; cf. Figure 1. The function 
id is an example of a member of *$>d that is supported on a spherical cap of radius 2r. If 
we would like to have a convolution root that is supported on a spherical cap of radius r, 
such as f,i(r) -1 / 2 l{#(-, •) < r} for id, it may not be suitable to choose all coefficients of 
the convolution root nonnegative. In the Euclidean case, the existence of convolution 
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roots with half-support, so-called Boas-Kac roots, is discussed in Ehm et al. (2004) 
building on the classical result of Boas and Kac (1945). It remains an open problem 
whether Boas-Kac roots always exist for functions in ^d- 



4 Differentiability 

4.1 Proof of Theorem 1.2 

We denote by ^>d the space of all continuous functions (p : [0, tt] — > R which are such 
that the function tp(9(-, •)) : S d x S d — > K is positive definite. The difference between 
the spaces ^d and *&d is that the members ip G ^d C *&d are additionally required 
to fulfil ip(0) = 1. Theorems 3.1 and 3.2 also hold for the class ^ with the obvious 
modification that we need to require X^nLo bd,n < 00 instead of X^nLo ^cf,n = 1 f° r the 
Schoenberg coefficients in the former. 

For the proof of Theorem 1.2 on the differentiability of positive definite functions 
on spheres we show the following proposition, which can be applied iteratively to yield 
the assertion. 

Proposition 4.1. Let d > 1, ip G ^d+2- Then ip is continuously differentiable in (0, tt) 
and its derivative can be written as 

W) = -^(/iW-/ 2 W), 

where /i,/ 2 £^ 

Proof. By DLMF, 18.9.19 the derivative of for a > and n > 1 is given by 

± c C«(x) = 2aC^\{x). (13) 
We assume first that d > 2. As $?d ^> we can write ip as 

n=0 l-U 

with non-negative coefficients 6rf in such that X^nio ^ rf > n ^ °°' see Theorem 3.1. For 
N G N, d G [0, tt] we define 



n=0 °" \ L ) 

By Proposition 2.7 ^jv converges uniformly to ip. Let § G (0, 7r). By (13), the derivative 
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of ipM is given by 



N r (d+l)/2 ( 

^ w = e ^(d - 1) (- sin ^> 

n=l 

_J_f^ 1 / (n + d-2)(n + d-l) (<i .. 1)/2 

n(n + l)_^ 



2ra+d-l n+1 

J_f, n(» + d - 1) / gfc^W) Ca i)/2 (cos^) \ 
"sin*^ d ' n 2n + d-l ^ ^ + - 1)/2 (1) C£i 1)/2 (1) / 

where we used (5), and 

(n + l)(n + 2) ^ + -i)/ 2(co ^ ); 



(d- l)(2n + d+l) 
see DLMF, equation (18.9.8). Therefore 

d 



(sin 0)^(0) = -b d . 



d+l 

" ( n (n + d-l)^ {n + 2)(n + d+l) u \ C r ( ^ 1)/2 (cos I) 
+ ^ V 2n + d - 1 d ' n 2n + d + 3 d ' n+2 J r (d_1)/2 m 



(n + 2)(w + d+l) C&i 

r (d-l)/2, , ■ 
n=N-l " ' ' °n+l 

The last term in the above equation converges to zero uniformly in i? as N — > oo by 
Gneiting (2012, Corollary 4) and Lemma 4.2. We will omit it in the sequel. Using 
Gneiting (2012, Corollary 3(b)), we obtain 



n(n + d-l) (n + 2)(ra + d+l) 

~°d,n ^ ; — T~ ~ "d,n+2 



2n + d-l " 2n + d + 3 



dn d(2n + d+ l)(n + 2) 

Od+2,n 77; , j , »w — 7T"d,n+2- 



n + d ^ ' (2n + d + 3){n + d) 
Hence, 

(smtf)^) = d^ 6^—^—— 

n=0 U n+1 \ l > 

^ (2n + d-3)n cg^^cos g) 

£j (2n + d-l)(n + d-2) d ' n ^_" 1)/2 (1) 
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We set $ ] = 0, 



) = d n , 1 : ferf+2,n-l, for n > 1, 

n + a — 1 



and 

_ f2) (2n + d-l)(n + l) , , 

(2n + a + l)[n + a — 1) 

The sequences {/3n^}neN > * = 1)2, are nonnegative and summable by assumption. 
Therefore they are the Schoenberg coefficients of some functions f\, f 2 E ^rf. By 
Proposition 2.7 their partial Gegenbauer sums converge uniformly, which yields the 
claim. 

If d = 1, the proof uses the same arguments with Gneiting (2012, Corollary 3(a)) 
instead of Gneiting (2012, Corollary 3(b)). The Schoenberg coefficients of the functions 
/i, f 2 are then given by fffl = {{n- l)/ra)6 3)ri _ 1 , /3i 2) = &i,„+i, for n > 1, and f3^ = 0, 
/3j 2) = (1/2)6!,!. □ 

Lemma 4.2. Let (a n ) ng N an increasing sequence converging to 1, s-uc/t i/iai the se- 
quence (a^) nS N «s bounded away from 0. Suppose that X^^Li b n < oo for some sequence 
(6 n )neN °/ nonnegative numbers. If 

b n > a n b n+ i, for all n G N, 

i/ien n b n — > as n oo. 

Proof. Let (a") n€ N be bounded below by C > 0. Let e > 0, choose no, such that 
Ylk=n+l bk < £ for all m > n > no- With m = 2n we obtain 

2n 2n 2n-l 2n 

S > b ^ E II ^ > E («n) 2n "^2n 

fe=n+l fc=n+l j'=fc fc=n+l 



> al n nb 2n > C 2 n6 2n > 0. 



Using the same argument for m = 2n + 1 yields the claim. □ 



4.2 Optimality of Theorem 1.2 

In this section we show that Theorem 1.2 is optimal for all odd dimensions using similar 
ideas as in Gneiting (1999). We are not aware of a function ip £ ^2 with discontinuous 
derivative. If such a function was available, our method immediately also yields the 
optimality of the differentiability result in even dimensions. 

We introduce a turning bands operator for isotropic positive definite functions on 
spheres in analogy to the Euclidean case, where the turning bands operator originates 
in the work of Matheron (1972). Let (3 = {(3 n )n£No be a sequence of real numbers. For 
an integer fceZwe define the sequence (3 o as follows. If k > its members are 

jo, i£n<k, 
(P ° n)n = < R 

[Pn-k, if n > ft 
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for n G No- If k < we put (/3 o Tk) n = /3 n -fc for all n G No- Let d > 1 be an integer. 
For a summable sequence (3 = (/3 n )neN of nonnegative numbers f3 n we define i^diP-,^) 
for $ G [0, 7r] as 

n=0 I 1 ) 

Proposition 4.3. Lei d > 1 be an integer and let (3 = (/3 n )neN ^ e a summable sequence 
of nonnegative numbers f3 n . Then, for all r G [0, ir], 

^d{/3,r) = /3 + cosr ip d+2 {P o r_i,r) + ^sinr ^ +2 (^ ° r), (14) 



and 



^(smr) d ^(^T.i.r) = jT (sin tf)^ 1 (MM) - fa)d&. (15) 



Proof. Suppose first, that d > 2. Using Proposition 2.7, (12), and (5) we obtain 



[ T (smV^MMW = f> f r (^) d - l Cn ^ 1) ( )^ ) d^ 
^° n=0 ^° C« (1) 

„ r 1 oo ^-,(d+l)/2/ 

= /3b/ (sin tf)^dtf + \ (sin r) rf ^ &, 



which implies (15). Differentiating both sides of (15) with respect to r yields (14). The 
case d = 1 can be shown using the same arguments. □ 

The proof of Theorem 1.2 shows that the differentiability of a function ipd{Pi ') only 
depends on the nonnegativity and the asymptotic properties of the sequence (/3 n ) ne ^ . 
Therefore, for any k G Z the function ipdift ° T kr) 1S continuously differentiable if and 
only if the same holds true for ipd(P, •)■ Let c G (0, it). Then the function 



=max{fJ, (l- i?g[0,tt 



belongs to the class Vl/i as can be shown by elementary arguments. Its first derivative 
does not exist at the point •& = c. Let (3 = (/3 n )neN be the sequence of 1-dimensional 
Schoenberg coefficients of ip. Let d > 3 be an odd integer. By (15) and the above 
remark on Theorem 1.2, the function ipd(/3 t_m-i)/2^) £ and its derivative of 
order (d — l)/2 does not exist at i? = c. 

The truncated power functions = max{0, (1 — $/c) T } were studied in detail by 
Beatson et al. (2011). They were able to show that they belong to ^ if r > (d + l)/2 
for d = 3, 5, 7 and conjectured the result for all dimensions. Theorem 1.2 immediately 
shows the necessity of the condition for all odd dimensions. 
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